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Abstract. In this article, the quantum Hamilton- Jacobi theory based on the position dependent mass model is 
studied. Two effective mass functions having different singularity structures are used to examine the Morse and Poschl- 
Teller potentials. The residue method is used to obtain the solutions of the quantum effective mass- Hamilton Jacobi 
equation. Further, it is shown that the eigenstates of the generalized non-Hermitian Swanson Hamiltonian for Morse 
and Poschl- Teller potentials can be obtained by using the Riccati equation without solving a differential equation. 
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1. Introduction 

Quantum Hamilton- Jacobi (QHJ) theory has always attracted much attention in both non-relativistic and relativistic 
quantum mechanics. An interesting work on the quantization of bosonic chiral Schwinger model is studied within HJ 
context pQ . QHJ theory, which is a generalization of classical Hamilton Jacobi formulation, based on transformation 
theory [2 . A fundamental work for obtaining the eigenfunctions of a quantum system with using the QHJ theory may 
be found in [3J. Furthermore, there has been an increased research interest in quantum trajectories in the complex 
plane [H[S]. The study of the De Broglie-Bohm approach is formulated in a time dependent domain [BJ. The QHJ 
formalism, first developed by Leacock and Padgett in 1983 [7]. Exact and quasi-exact solutions of a group of soluble 
potentials are also studied [8] as well as the supersymmetric quantum mechanics and shape invariance approach [5]. 
Moreover, Sanz et al showed that the chemical reactions are analyzed within generalized HJ framework |10j . Some 
detailed and interesting work on extensions of analytical mechanics to the complex plane can be found in . In [12] , 
the spectrum for the class of potentials within QHJ is obtained. 

On the other hand, position dependent mass model has been a subject of considerable interest in recent years 
[13TI22"] following the work of Von Roos and Levy-Leblond |TT1[20] that claims to describe some physical phenomena. 
The importance of this kind of Hamiltonian arises from the motion of electrons and holes in semiconductors that can 
be described by them. In this model, the mass is taken as the effective mass of the particle, which depends on the 
material |17) . Von Roos showed that Hamiltonians with position dependent mass are not Galilean invariant within the 
extension of Bargmann's theorem |20j . Later, Levy-Leblond showed that Galilean invariance can be used to obtain 
Hamiltonian of an electron with spatially dependent mass, a model of a hetero junction is proposed in |17j . 

Recently, a great attention has been devoted to investigate the properties of non-Hermitian Hamiltonians in 
theoretical physics [53JH1]. A wide class of VT symmetric, non-Hermitian Hamiltonians provides entirely real 
spectra [55] . For recent developments please see [26ff2"9"j . Mostafazadeh has investigated the larger class of Hamiltonian 
models known as pseudo-Hermitian, defined by the relation H' = ijHrj -1 , for a Hermitian operator rj [301132] . Various 
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papers appeared on this topic |33H36] . As a non-Hermitian and "PT-symmetric Hamiltonian model, some properties 
of Swanson Hamiltonian are studied in the literature 23,28,33,36 . 

In the present paper, as a first study of the QHJ equation introduced in a mass dependent form that can be 
a generalization of QHJ equation, exponential and hyperbolic type effective mass models are considered and the 
transformed equation, which is a Riccati type, is mass dependent. The effective mass-QHJ equation is solved for the 
effective potentials using a Laurent series. Moreover, a pseudo-Hermitian Hamiltonian model is connected with QHJ 
equation and two types of effective potentials are generated. 
2. Preliminary Exploration of Two Approaches 

In this section we shall introduce a brief summary of the QHJ theory and the effective mass problem. 
2.1. Quantum Hamilton Jacobi theory 

In [7J, the authors considered a particle moving in one dimension under the influence of a potential U(x). Classical 
transformation equations p — 4^, Q — are used in the Schrodinger representation which is H = H(x,p) = 
p 2 + U(x) where x, p are linear coordinate and momentum operators. Then, the QHJ equation is stated [7J as 



where W(x, E) is the generating function. If the quantum momentum function p is defined by p(x, E) = dW(x, E)/dx, 
then the QHJ equation is defined by 



There can be a transition between the Schrodinger and the QHJ equations via the wave-function defined by 



two solutions. But we have to look for a physically acceptable solution. The appropriate boundary condition was 
proposed as p(x,E) — >h^o Pci(x,E) in Leacock and Padgett's paper; this is only for simple potentials, because 




p 2 - ihp - (E- U{x)) = 0. 



(1) 



= en 



IW(x.e) jjj^ wn i cn satisfies the physical boundary conditions. ([T]) is a nonlinear equation and it has 
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Pd(x, E) = \J E — U(x) may have several branch points. The quantum action variable is defined [7] by 

J(E) = — I p(x,E)dx = rih, (2) 
2tt Jc 

which connects eigenvalue J to the energy E. Here, C is the closed contour which encloses the cut of p c i(x, E) between 
Xi and x%. Poles of p{x,E) are enclosed by the contour C, thus J{E) = nfi [7|. As it is seen from quantum action 
variable and boundary condition that p(x, E) has poles of residue —ih on the real axis. The details of the model and 
calculations for obtaining the spectrum will not be given in this paper, however the reader can look at [HOE] for more 
details. 

2.2. Effective mass problem 

Let us start with the general Hermitian effective mass Hamiltonian which was suggested first by von Roos [20 

H = ^(m a (x)prn^ (x)pm' y (x) + (x)pm 13 (x)pm a (xj) + V(x), (3) 

with the constraint a + j3 + 7 = —1 on the parameters, and m(x) is the position dependent effective mass. One can 
express this Hamiltonian in various ways, depending on the choice of the parameters. Thus, various special cases of 
the kinetic energy operator can be given |15j by : 

r = Jf- P 2 + P 21 ) (4) 

4 \m m J 

H( p s") (5) 
t-U-Lt-L). (e) 



2 \ yrn a/W/ 

In this work, we will keep the same general form of the Hamiltonian. Thus, the one-dimensional effective PDM 
Hamiltonian is given [20] by 

1 d2 M '( x ) d „ / x 

where 

1 M" a a AM') 2 



ff (x) = V(x) + -{fi + i)—- [ a (a + (3 + 1) + (3 + (8) 
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a, f3 are parameters as written before, and primes stand for the derivatives with respect to x. Note that we use the 

dimcnsionlcss form M(x) of m(x) = moM(x) and Ti — 2mo = 1. Now we can introduce the corresponding eigenvalue 

equation for ([7]) as 

1 ip"(x) M\x) ip'{x) 
M(x) ip{x) M z {x) ip{x) 

where e and (p(x) are the eigenvalues and eigenfunctions of the Hamiltonian ([3]). Our task is now to adapt © to QJH 
theory. 

3. Quantum Effective Mass- Hamilton Jacobi Model 

Now we make a new transformation as p(x, e) — — i^— in (|9|). Hence we have, 

, m'(x) 
-ip (x,e)+p 2 (x,e)+i—^p(x,e)+M(x)[V eff (x)-e}=0 (10) 

that will be called as the quantum effective mass(QEM)-Hamilton Jacobi(HJ) equation. We shall use p for quantum 
effective mass momentum function. It is obvious that the QEM-HJ equation can be transformed into the QHJ equation 
by using 



where <p(x) — y M(x)ip(x). And we also point out that tp corresponds to the eigenfunctions of the Schrodinger 
equation. Thus, the relation between two quantum momentum functions is 

M'(x) , . 

Introducing a general form for p: 

p = a(x) + b(x)p + c(x)p 2 , (13) 

which has the same form as (fT0|). and using 

-m 1 i \ 6 ( x ) c '( x ) nA\ 
V{x) = -r^v(x) - —— - 14 
c(x) 2c(x) 2c^(x) 
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in (fT5|) . we have then, 



v (x) -v 2 {x) =G(x), (15) 
where G{x) is given by 

b 2 b' be 3 c' 2 c" 

^)=-- T +2-2^-4^ + 2c' (16) 

and primes denote derivatives with respect to x. If we compare (|15[) . (| 16[) . (|13l) with (|10p. we can express a, &, c in 
terms of the physical quantities: 

a(i)=-iM(.)(K / /(x)- £ ) (17) 
Af'(x) 



b(x 



M(x) 



Then, G{x) can be written as 



4M 3 (x)(V efJ ix) - e) + 3M' 2 (x) - 2M(x)M"(x) 
4M 2 (x) 



j4 short discussion on the multiplicities 

G(x) may have the number of poles because of the mass and effective potentials. In the previous works, there were 
not any comments on the multiplicity of the poles [3,8,9, 12 . In terms of the multiplicities of the poles of G(x), the 
maximum number of distinct meromorphic solutions can be discussed. 

Theorem : If G(x) has at least one simple pole, then (1151) admits at most one meromorphic solution. 

Corollary : Let G(x) is a meromorphic function, then it has at least one pole. If (|15[) has two distinct meromorphic 

solutions, then all poles of G{x) are of even multiplicity. 

Proof of Theorem: Let v be a meromorphic solution of (|15|) and xo is a simple pole of G{x). Obviously, the residue of 
v at Xq is — 1. Thus, there must be a neighborhood of xo, which is T>, such that the function 

£ := {x-x )e~f*o xdx (19) 
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and the function x are analytic in a neighborhood of Xq, and \ satisfies \ = v + x -x ' Henceforth, £ satisfies the 
well-known differential equation 

C" + G(x)£ = 0. (20) 

Suppose that the theorem is false and (|15p has two distinct meromorphic solutions Vi, Vi- Let £i, £2 and \\i X2 satisfy 

(f20f and (fT9|) on a neighborhood 5. If £1 and £2 are linearly independent, then; ^ (£1757 — £2^) = 0. From this 

point of view, there is a function / and a constant c such that / := f 2 -, / := then xo is a simple pole of /. By 

si ?i 

using (|19l) . / becomes analytic at xq. And there is a contradiction. □ 

Here, Theorem 6.1 in Ref. |37j is followed to introduce the Theorem given above. In fact, previous works include 
even multiplicity. But, there may be also simple poles depending on the mass function. We will see this in the example 
of the Poschl- Teller potential, discussed below. 

3.1.1. Singular Mass with simple pole: xq — ikir, k — 0,1,2,3,.. We now discuss the QEM-HJ model for 
M(x) = ■ \ . We introduce 

V(x) = Ficothx + F 2 + F 3 sinha; (21) 

then we obtain G(x) and V e f / as 

- G(x) = Vi cothxcscfcr - (a(a + /3 + 1) + -) esc h 2 x + - +V 3 + ^— ^ + - (22) 

4 sinhx 2 4 

V e ff(x) = Vi cothir + V 2 + [V 3 + ^-y^ + ((3 + l)csch 2 x- (a(a + f3 + 1) + + 1) coth 2 x] sinhx. 

Thus, a Riccati equation can be obtained by the appropriate transformations (given above) as 

v - v 2 - V\ coth x esc hx + {otiot + /3 + 1) + -)csch 2 x - 2 ~ £ - V 3 - ^— i -1 = 0. (23) 

4 sum 2; 2 4 

There is a simple pole coming from the coefficient of V% — e in the above equation. This equation is solvable for V% = e 
as a special case. On the other hand, one can expand G{x) in a Laurent series as 

G(x) = ^±^ +*»+... (24) 
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where k% = —V\. k% = (a (a + (5 + 1) + i and /c 3 = e — V^. There is one simple pole and also even multiplicity. We use 

the mappings, v = — \Jy 2 — \ (c(y) — 2(y%-i) ) > V = cosn:c > in (|23|) and we find 

j , , 2 _ 1 , 3t/ 2 1 / Fiy _ q(a + /3 + l) 1 + 1 

^ ^ 2(j/ 2 -l) 4(y2-l)2 + (y2_l) i v (y2_ 1) ( y 2 _ j) + "3 f 4 + 2 



Now we have to seek solutions to the above equation. The first task is to expand £ in Laurent series and to find the 
residues. Then, £ is expanded in a Laurent series as 

C=-^r+Po+Pi(tf-l) + ... (26) 

and the residue b\ at y = 1 can be obtained, 

1 



6i |„=i = - (l ± Vyi-a(a + /3 + l) + lJ . (27) 
At y = — 1, the residue is 

bi\ y =-i = i (l ± v 7 -^! - "(a + /3 + 1) + l) ■ (28) 

Because we cancel the term fc 3 , there are two solutions for each value of the residue, as mentioned in the Theorem. 
But, because of the integrability condition, we will choose one of them. Now, F (y) = Q(y)F(y) will be used in 

y-1 y+1 

where h is a constant because of Liouville's Theorem [38 . This equation will be put in (|25[) to obtain the spectrum. 
The constant h — can be obtained here. In the limit of y oo, Q(y) behaves as Q(y) — > *J and Q + Q 2 — > "^i" 1 ^ . 
From the coefficient of y -2 one obtains 

2&i| 1/ =i6i| 1/ =_i + 2n(6i|j,=_i + 61^=1) + n(n - 1) + | - V3 - a(a + + 1) = 0. (30) 

o 

Hence, from (1501) we have 



V 3 = 7 2 ~ + V-^i " "(« + P + 1) + 1 + VVi ~ "(« + P + 1) + 1) - «) 2 , (31) 

where 

+ 1 1 1 



7=-^ + -A-l + V-Vi -a(a + /3 + l) + l + ^Vi-a(a + /3+l) + l) 2 . (32) 



2 4 4 
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We note that the negative sign is used for the residues because of the integrability condition. 
Finally we can give the solutions by using (l27l) . (f28|) . (|26|) and <p = e l 'J 



lMtf) = (l/-l) Bl -* (y + lf 2 -^ ^ § ' 2Q2_I) (2/), (33) 

and 

Mi/) = (i/-l) ai ~ s (y + l) 02 --* P^ 1 ' 1 ' ' 2a2 '~ 2 \y), (34) 

where P„'(y) are Jacobi polynomials and the constants are given by 

ax = |(1 - 2y/-V! - a(a + [3 + 1) + 1 + y/Vx - a(a + (3 + 1) + 1) 

a 2 = + - a(a + /3 + 1) + 1 + V^i - "(« + /3 + 1) + 1). (35) 

Poschl- Teller potential is generally given by [31] 

= A 2 + ( J B 2 + A 2 + A)csc/i 2 x- J B(2A+l)cothxcsc/ix. (36) 
Thus, if we compare G(x) in (|2"2"| with we can obtain the constants A and -B in terms of the constants of (|2"TT) 

A = - (-1 + ypV\ -a(a + /3 + l) + l + ^/Vi - a(a + /3 + 1) + l) (37) 

S = - (V^i -a(a + /3 + 1) + 1 - ^/-Vi -a(a + /3+ 1) + l) . 

3.1.2. Exponential mass with the pole: xq = — oo: We will use effective mass M(x) = e~ 2x in this case. In |13j . the 
authors introduced V(x) as 

V(x) = V e 2x - B(2A + l)e x (38) 

and G(x) and V e ff have the forms 

G(x) = ee~ 2x + B{2A + l) e - x -V + 2(/3 + 1) + Aa(a + + 1) + 1 (39) 

V eff = (V - 2(13 + 1) - Aa{a + (3 + l))e 2x - B{2A + l)e x . 
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If we use the mapping p — i(l + v), v — y(£ — i) given in (TT4"]) and (IT71) and y = e _:c , we have 

C + C 2 + + ^ (e?/ 2 + S(2A + l)y - v + C) = 0, (40) 
where C = 2(/3 + 1) + Aa(a + + 1) — 1. Now we expand £ in a Laurent series as 

C= —+PQ+P1V + - (41) 
V 



and the residue 6i became 6! = i(l ± 2-y/| Vo — C\). We now discuss the residue at y — oo, on the extended complex 
plane. Using y = j in (|40|) and expanding £ = B a + Bit + B 2 t 2 + ... one obtains the residue i?i = — B ^i?o^ ' 
-Bo = —\/—e. The coefficient of y -1 for large y gives 

U = 2(/3 + 1) + 4a(a + /3 + l)-l + (A- n) 2 , (42) 

where £ = — B 2 [13]. The solutions can be obtained following the same way as it is discussed before, if C is written 
C = y + Q(y) + h, and F (y) — Q(y)F(y) is used, then the constant h — —\ can be obtained. Using (|4"0"]) and the 
residue b\ one obtains 

Vn(y)=y A - n+1 e-^L 2 j A - n Hy), (43) 
which is consistent with [13] . 

4. Pseudo-Hermitian Hamiltonian connected with QEM-HJ model 

Suppose that H is not a Hcrmitian Hamiltonian. Then, 77-pseudo-Hcrmiticity of H is equivalent to the condition given 
by [3DJ 

i? 1 = r/iJr/" 1 , (44) 
where 77 is a linear, invertible operator. The metric operator 77 is not unique for a pseudo-Hermitian operator H. 
There exists a mapping from the non-Hcrmitian H to its Hermitian counterpart h, through a similarity transformation 
h = pHp -1 , here H, a differential operator act in a complex function space and the mapping function is positive definite 
p = y/rj |30] . A non- Hermitian, VT symmetric Hamiltonian model is considered by Swanson and it is given [36] as 

H = cj(a t a+i) + aa 2 +^a t2 , (45) 
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where a, generalized creation and annihilation operators, 

. , s d , , + , , „ d , , oVl(x) 

a = A(x) — + £(x), a f = -A(x)— + B(x) 

ax ax ax 

and cj, a, /3 are real, positive constants. In [24], the more general form of Swanson Hamiltonian is introduced as 

H = w(a t a + i) + aa 2 + 13a) 2 + ja + Sa\ (46) 
which can be written in a differential operator form |24j 

where u> — u — a — /?, 

ftifar) = (a - /3)A(x)(2B(x) - A (x)) + ( 7 - <5)A(x), (48) 

and 

c 2 (x) = (w + a + /3)B 2 (x) - (w + 2/3)A'(x)B(x) - (w - a + /3)A(x)B' (x) (49) 
+ /3(A(x)A" (x) + A' 2 (x)) + ( 7 + <5)B(x) - <L4'(x) + ~ 
Let '5 and £ are the eigenfunctions and the eigenvalues of (|4"Tf i.e., 

= (50) 

If p is introduced by 

and if it is used in (|50|) . then it yields 

p (x) = -i P 2 ( X ) + hix) ^f {x) p( X ) ~ ■ (52) 

Comparing above equation with (|13j). we obtain 



&(x) 



C2 — £ 
' QA(x) 2 
bi - 2ujA(x)A' (x) 



a(x) = — i — 2 (53) 



wA(x) 2 
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Meanwhile, it is obvious that (|5"2"j) is not any type of the QHJ equation that is discussed in Sections 1 and 2, because 

(|47|) is not Hermitian, and H 1 ^ = £^ is not a regular Sturm-Liouville type differential equation. But one can find a 

Hermitian counterpart of (|47l) by using a convenient metric |24j . But here, we can find a Riccati- type equation which 

corresponds to a QEM-HJ equation. We then use: 

in order to transform (j5"2")) into (TT5")) . Thus, G(x) can be obtained as 

= A"(x) b[(x) _ c 2 {x)-£ _ bj(x) 
[} A(x) 2u>A 2 (x) ljA 2 (x) 4w 2 A 4 (j:)' 1 ' 

We note again that, the primes in equations given above denote derivatives with respect to x. Hence, we can introduce 
^>(x) as 

^(x) = e J '"(''I' (f>(x) (56) 
where 4>{x) = e V ) J _ 

On the other hand, if we follow a standard procedure(for the Schrodinger operator) [24,34,35 , f|50[) can be 
transformed into an equivalent Hermitian form by means of the similarity transformation |30) 

h = pHp~ x , (57) 
where the mapping function reads [24ll34l l35] 

p = exp(-± / dx^r) (58) 



2lo J A 2 (x) 
such that rj = p 2 . This transformation leads to 



where [24] 



h = ~ Q iL A2{x) iL +Veff{x) (59) 



V ef/ = ( ( " _ P)2 +uj + 2{a + p)) B(x)(B(x) - A (x)) - (u> + a + f3)A(x)B(x)' + ^±Ia(x)A(x)" 



w + a + /3 



(60) 
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4-1. Poschl-Teller and Morse potentials in Swanson model 

Now we will study the effective potentials in Swanson model. If we compare |59|) and ((7]) it is seen that 



Let us discuss two different types of A(x) function. If we can give an ansatz for A(x) and B(x), we may be able to 
obtain some effective potential models. In this manner, one can obtain the residues and all solutions in terms of the 
Swanson parameters uj, a, ft. 

i) Firstly, we take 7 — S — for convenience. If we choose 



A(x) = J^, B(x) = A(xy x 



then, using ([551) . we obtain 



G(x) = - " 2 + ^-^?-MP + ^) _ ^28 ^ + (62) 



cothxcsc/ix - f^- + w - 4a/3 csch 2 x. 



2oj 2 V 16w 2 4 

and from ([SU|) we have 

T4/ / (x) = — + ^1 jl cot _|_ _^ — gQ^jj x csc ^ x _|_ ^2 _ csc ^ x (63) 
2w 16a; 2 

Q + • , 
4- — ; smhx. 

4w 

The QEM-HJ equation is solvable for u = 2£. Comparing (|62"|) , (fl~5)) and ([36]) we get 



(a + /3)(q + /3-3^)+^ 2 + 4a/? 

= 2a(a + ^-o,)» (64) 

(a + /3) 2 - w(a + /3)(a + 3) + u? + 4(a - l)a/3 

£> = — 



2a(a + /3-w) 2 

Vi, V2 and V3, the parameters of (|38|). can also be written in terms of Swanson parameters as 



Vi= _ ^a+f)-4al 
V 2 = w 

a 2 + (/? - 2w) 2 - 2a(j9 + 2w) 1 /3 + 1 



^3 = 



16w 2 4 2 
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e =2£, 

where j3 is used instead of j3 that is the parameter given in von Roos Hamiltonian. The solutions ^ which is defined 
in (|56[) is then given by 

~ (coshx - 1) ^ so (coshx + 1) 2 sa 2 - 2; (coshx), (66) 

where the constants A and B were given by and < x < 00. In this case we note that the metric is given by 

1] = (sinha;) - ^" ^tanh — J (67) 

such that the Hermitian inner product <C ^f(x)\^(x) ^> r) =< $(x)\r)\ i &(x) > can be used for the normalization. The 
constraints for the parameters are given by B — A > ^'^g^ 4 "^ which is satisfied by the relations A and B in (|64[) . 
ii) Next we introduce A(x) and B(x) as 

A(x) := = : B ( x )- 

In this case we obtain G{x) as 

rw . (a-/?) 2 2(/3 7 + ^)-( ( 5 + 7 y / w ( 7 -^) 2 \ -2x 



and the effective potential is given by 



By comparing and we obtain 



and 



A= _l_ + 2W, + aS)-(S + 7 ) B _, 
2 4w5 

CJ ( 7 -,5)2\2 

B -'^2 +1I 4^ 1 
e = £ 



= 2i . ; 1 1 + 4,i . n - ; 1 1 - :i ( SLJ^j * , = A - n, (71) 
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where a, /3 are the constants defined in the von Roos Hamiltonian instead of a, j3. And the complete solutions are 

given by 

* n ~ {e- x ) A - n -^ e ( - B+ ^5 )e ~* L^ A - n \2Be~ x ). (72) 

In this case, r\ can be given as 

V = e -^ x+ ^ e (73) 

which can be used in the Hcrmitian inner product <C ^(x)\^(x) 3>^=< ty(x) \rj\^(x) > to obtain normalized wave 
functions. We note that the boundary conditions imply A > n and — oo < x < oo. 
5. Conclusions 

In this work we have studied the generalized QHJ equation within the context of effective mass approach. We have 
used p instead of p as an effective mass-quantum momentum function which can also be used to obtain the solutions 
of a position dependent Schrodinger equation, without solving it. We have obtained the solutions of the QEM-HJ 
equation for two different types of mass function. It is seen that if there is an even multiplicity of the poles, there are 
at most two solutions but only one of them is acceptable by reason of the physical conditions. On the other hand, 
the mass dependent QHJ equation differs from the regular QHJ equation because of the pole contribution of the mass 
term. This can be observed in the Poschl- Teller example such that the coefficient of the single pole, k^, is taken as 
zero. But, in the example of the exponential mass, the pole structure of G{x) does not lead to such a restriction. 

We have emphasized that the transformation in (|14|) that maps the QEM-HJ equation into a normal Riccati 
equation is also used for the Swanson model. This also means that an eigenvalue equation of the non-Hermitian 
system can be mapped into a Hermitian one by (|14[) which is seen to be in Riccati equation form. We have also 
shown that, using certain choices of A(x) and B{x), leads to some effective potentials generated from the Hermitian 
counterpart of the Swanson Hamiltonian. The solutions of (|52[) are obtained for the Poschl- Teller and Morse potential 
cases without solving the QEM-HJ equation. To obtain the normalized wave functions, we use the r\- Hermiticity 
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condition for (|6"6")l and (|72p . On the other hand, it is clear that an extensive class of soluble potentials can be classified 



according to choices of A(x), B(x). It would be an interesting problem to discuss a relativistic quantum stationary 
Hamilton- Jacobi equation within effective mass approach. 
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